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(Dated: February 1, 2008) 

The effect of strong shear flow on highly fluctuating lamellar systems stabilized by intermembrane 
collisions via the Helfrich interaction is studied. Advection enters the microscopic equation of motion 
for a single membrane via a non-linear coupling. Upon coarse-graining the theory for a single bilayer 
up to the length scale of the collision length, at which a hydrodynamic description applies, an 
additional dynamical coupling is generated which is of the form of a wavevector-dependent tension 
that is non-linear in the applied shear rate. This new term has consequences for the effects of strong 
flow on the stability and dynamics of lamellar surfactant phases. 

PACS numbers: 61.30.Dk Continuum models and theories of liquid crystal structure, 05.10.Cc Renormaliza- 
tion group methods, 87.16.Dg Membranes, bilayers, and vesicles 



I. INTRODUCTION AND OVERVIEW 

Dilute solutions of lamellar phases typically consist of 
highly fluctuating layers. The wide equilibrium layer 
spacings are governed by the interplay between the long- 
range steric repulsion, known as the Helfrich interaction 
and the bending elasticity of the bilayers. When these 
systems undergo flow, a range of interesting phenomena 
is observed, transitions to multilamellar vesicles (|, |j| 
and a reduction in layer spacing Q. Unlike layered one 
component melts, such as thermotropic smectics or di- 
block copolymers, flow can have a significant effect on 
the microstructure of the layers. Although flow certainly 
stretches the chains in diblock copolymers Jp| and can 
induce layer tilt in thermotropic smectics fl, the effect 
on the highly-fluctuating many-component layered sur- 
factant phases should be much more dramatic. 

As an initial step to account for some of this flow be- 
haviour, we previously conjectured that flow induces 
an effective anisotropic tension parallel to the flow, when 
lamellae are aligned in the c orientation (Fig. |l|). The 
"tension" is a response to projected area changes and 
acts to suppress the fluctuations. This led to predictions 
for either changes in layer spacing or an undulation in- 
stability. In a related work, Zilman and Granek H also 
proposed an effective tension, but isotropic, negative in 
sign, and of a different physical origin. While both stud- 
ies relied on inserting the "tension" heuristically into the 
dynamics as an effective free energy term, it is of interest 
to examine the dynamics of a membrane in flow to see 
how such a response can be generated dynamically. In 
this work we consider a lamellar phase stabilized by the 
Helfrich interaction, neglecting the effect of electrostatic 
forces. By coarse-graining the dynamics up to a length 
scale characteristic of the long wavelength hydrodynamic 
description, the typical transverse length L p between col- 



lisions, we demonstrate that flow can indeed induce a 
dynamical suppression of fluctuations that resembles a 
wavevector-dependent "tension" . 
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FIG. 1: The allowable steady state orientations a and c of a 
lamellar phase in uniform shear flow. 



The linearized relaxational dynamics of a membrane 
parametrized by a single height variable h subject to an 
anisotropic tension a, written in Fourier space, is 



d t h(q) = -A(q) + aql] h(q), 



(1.1) 



where K is the bending modulus, q is the wave vector and 
the kinetic coefficient A(q) depends on the particular re- 
laxation mechanism. The tension penalizes fluctuations 
and, if applied to a stack of such membranes that in- 
teract via collisions, would change the "preferred" layer 
spacing and renormalize the coarse-grained smectic layer 
compression modulus B. Our task here is to derive an 
equivalent term that would contribute to the effective dy- 
namics of highly fluctuating membranes in shear flow in 
the c orientation, v — jzx. We shall find that the equa- 
tion of motion for the coarse-grained height field in shear 
flow becomes 
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d t h(q) + ii^2(q x - k x )h(k)h(q - fe) = - {A{q)nq 4 + A x q 2 x ) h{q), (q > 



(1.2) 



r 



where the function A x depends on the wave vector, the 
strain rate 7, the kinetic coefficient A(q), and the coarse- 
graining length. The non-linear advection term, which 
arises from assuming an affine distortion of the mem- 
brane, induces the non-trivial renormalization of the dy- 
namics under coarse-graining. Hence, the "tension" , as 
such, is given by 



a = A x A- 1 (q) 



(1.3) 



and is a function of the strain rate and is typically wave 
vector dependent. 

Our calculation is similar to renormalization group cal- 
culations of the anisotropic Burgers equation, studied in 
the context of sandpilcsS, and of phase separating sys- 
tems under shear flow ]10|, |11| . However, unlike the calcu- 
lation of critical behavior, in which coarse-graining pro- 
ceeds until scaling is found, coarse-graining in our case 
can only be performed up to the natural physical cut- 
off corresponding to the collision length L p . This may 
or may not be in the scaling regime; if it is not, then 
there are generally other terms generated by the coarse- 
graining procedure. However, such a calculation is gener- 
ally impossible, so instead we estimate the coarse-grained 
dynamics from that obtained by a non-trivial scaling 
fixed point calculation. Moreover, we expect that terms 
that eventually become irrelevant in the scaling regime 
are in the process of being driven close to zero at the true 
coarse- graining length L p , and in any case are of higher 
wave vector and hence not present in a hydrodynamic 
description. 

Our one loop perturbation expansion of the dynamics 
yields, for any A(q), a tension that scales as a ~ j 2 . 
By considering the energetic cost involved in bending 
and stretching a single membrane, Zilman and Granek 
H estimated a tension (isotropic) with the same scaling. 
However, upon performing a coarse-graining one finds a 
dependence a ~ 7 e , where e depends on A(q). 

In Section [n| we summarize single membrane relax- 
ation dynamics, discuss previous studies of the dynamics 



of h and the coarse-grained smectic displacement vari- 
able u in shear flow and derive the equation of motion 
for a membrane in shear in the c orientation. In Sec- 
tion III we outline the coarse-graining procedure and ex- 



tract results for the renormalized dynamics, for different 
relaxation mechanisms. We conclude in Section IV with 
a discussion. The appendices collect calculations of the 
permeation dynamics and the details of the renormaliza- 
tion calculation. 



II. MEMBRANE DYNAMICS 



A. Membrane Dynamics without Shear Flow 

First we review the equilibrium properties of a 
Helfrich-stabilized lamellar phase, in preparation for 
studying its dynamics in flow. The long-ranged entropic 
interaction, characteristic of such a phase, is a conse- 
quence of colliding membranes due to thermal fluctu- 
ations. A key notion is the characteristic distance L p 
between the collisions, marking the transition in length 
scale between the membrane and bulk smectic behaviour. 
This enables us to relate the static behaviour of a single 
membrane at a mean layer spacing £ to the compression 
elasticity of the lamellar phase. 

In general, if h(x, t) is the height of a membrane above 
a plane in d dimensional space, where x — (x, x\) is a 
(d — 1) dimensional vector in the plane and t is time, the 
elastic free energy in the Monge gauge ({(Vh) 2 } <C 1) is 
given by the Helfrich Hamiltonian 



mgc 
1 



(2.1) 



where the bending 

\3-d 



modulus K has dimensions of 
(energy)*(length) a_d . For the rest of the paper we use 
discrete or continuous Fourier transforms as convenient, 
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(q-x — uit) 



(discrete) 



h(x,t) = < 



I 



(2.2) 



where q = g x ii + • x± and is the frequency. For 
clarity we reproduce results in the rest of this section for 
d = 3 so that xj_ = yy. In terms of the Fourier com- 
ponents h(q), the equipartition theorem gives the equi- 
librium height correlation function (\h(q)\ 2 ) = T/K,\q\ 4 , 
where we take the Boltzmann constant ks = 1- The 
statistics of small fluctuations with wave vectors qL p > 1 
are adequately determined by the Helfrich Hamiltonian. 
However, the behaviour of large fluctuations with wave 
vectors qL p < 1, is complicated by the steric repulsion. 
Constraining the fluctuations to within a layer spacing 
determines the mean transverse length L p between col- 
lisions, which scales as L p ~ I^Jk/T. The loss of en- 
tropy associated with each collision contributes to the 
free energy change of the lamellar stack under a change 
of layer spacing, and hence a bulk compression modu- 
lus that scales as B ~ T 2 / n£ 3 jjj. Long wavelength 
{qLp < 1) lamellar behaviour is described in terms of 
the displacement u of the mean layer position from its 
mean value, defined by 



u(x, y,z — n£) 



dfx' 

[h(x' -x,y -y)- n£] 
A{L P ) a V l p) 



.(2.3) 

where A(L p ) is the typical membrane area per collision. 
Bulk equilibrium smectic behavior can be obtained by 
removing the small wavelength, high q, height degrees of 
freedom h, until a coarse-grained description entirely in 
terms of the long wavelength variable u remains ]T^ | . Our 
goal is to carry out this procedure for the dynamics of a 
single fluctuating membrane in flow, by coarse-graining 
the dynamical description up to the collision length. A 
complete calculation would naturally need to simultane- 
ously incorporate the steric repulsion. 

The relaxation dynamics of a single membrane can be 
described by a Langevin equation, 

d t h(x,t) = - fd 2 x'A(x-x') J. H +$(a,t), (2.4) 
J dh(x',t) 

where the thermal noise £(q, t) describes the neglected 
microscopic degrees of freedom. The noise has zero mean 
and a variance given by the Fluctuation Dissipation The- 



(£(<7i, ti)Z(q 2 , ta)) = 2TA( Ql )6 ( Ql + q 2 ) 5 (h - t 2 ) 
A spatial Fourier transformation gives 

5H 



d t h(q,t) = -A(q) 



5h(-q) 



(2.5) 

£(q,t). (2.6) 



The kinetic function A (q) depends on the details of the 
fluid-membrane coupling; a general form is 



A(q)=A |g|' 



rf- 1 r n+1 \q\ m , 



(2.7) 



where the exponent m and the associated length scale I 
depend on the relaxation mechanism. 

Three relaxation mechanisms are summarized in Fig. |^ 
with relaxation functions given by 

ri~ 1 q~ 1 m = —1, isolated 
A(q) - <( r]^ 1 (q° m = 0, permeable (2.8) 
ri~ 1 l 3 q 2 m = 2, confined fluid. 

Thin, impermeable membranes exhibit two regimes: (i) 
(m = —1) For q^ 1 > £ the membrane is damped by vis- 
cous solvent drag. Solving the linearized Navier-Stokes 
equations for the solvent flow yields A(q) = l/Ar/q [ p^[ , 
whence iu> — nq 3 /4r). (ii) (m = 2) For g _1 < £, solvent 
flow is screened by the surrounding membranes, leading 
to iu> = K£ 3 q 6 /16r/ (w). Permeable membranes exhibit an 
additional regime, (iii) (to = 0) For £ < q^ 1 < £ where £ 
is a permeation length scale that depends on the size and 
the density of microscopic defects and membrane thick- 
ness, A(q) ~ e/77. A simple model of cylindrical pores 
(common in lamellar phases) of width w and mean sep- 
aration R leads to £ ~ R 2 t/w 4 , where r is the thickness 
of the membrane (see Appendix [X]). We may envisage 
particular systems in which, at sufficiently high perme- 
abilities < £~ 1 the isolated regime is excluded, and 
others for which < L p 1 so that only the permeable 
regime remains. 



J 



B. Equations of Motion in Shear Flow equations of motion previously used to describe lamcl- 

1. Previous studies: linear advection 



Before we consider the dynamics of a single fluctuat- 
ing membrane in flow in the c orientation, we review the 
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Confined Fluid 



Isolated 



Permeable 



i 1 r* C" 1 

A(q)= 4 2 /ri A(q)= l/r\q A(q)= 



e-i 



FIG. 2: Different scaling regimes for the kinetic coefficient, 
dynamical description breaks down. 



( is the permeation length and e is the length at which the 



lar phases in flow. Milner and Goulian |15| explored the 
stability with respect to a and c orientations of a ther- 
motropic smectic under flow. Since the layer fluctuations 
of most thermotropic systems are small relative to the 
length scale of the layer spacing, shear will not have a 
significant effect on the shape of the layers or the inter- 
nal structure. Thus the appropriate parametrization is 
entirely in terms of the broken symmetry smectic dis- 
placement variable u(r). In flow, this quantity simply 
advects as a passive scalar. In simple shear flow with 
average flow velocity parallel to x, and velocity gradient 
direction n, the flow field is 



v(r) = (r ■ n^x 

n = cos 8z + sin By. 



(2.9) 



where fi = y in the a orientation and n = z in the c 
orientation. In this case, the dynamics of the smectic 
phase is 



\ - "fq x n 



d_ 
dq 



u(q,t) = -13(g)- 



5F 



5u(q,t) 

(2.10) 

where T is the smectic free energy, \{x, t) is the noise 
and j3(q) is the smectic kinetic coefficient. After minimis- 
ing an effective free energy expanded in powers of strain 
rate, Milner and Goulian showed that the a orientation 
is stable with respect to the c orientation; moreover, they 
demonstrated that a more rigorous calculation of the dy- 
namic response function incorporating non-linear terms 
in the free energy exhibits the same result. 

Earlier, Bruinsma and Rabin [[l6) studied the effect of 
shear on a lyotropic smectic phase in the c orientation. In 
most of their calculations they assumed that flow doesn't 
change the membrane shape and character, and hence 
considered the same passive advection of u when calcu- 
lating the effect of shear on the smectic hydrodynamic 
dispersion relations. However, they did consider the de- 
gree to which flow influences the microscopic height vari- 
able h, and estimated the shear rate j c at which an affinc 
deformation of individual membranes leads to significant 
suppression of fluctuations, 



7c 



rpb/2 

K 3/2 r jg3 



(2.11) 



This was obtained using, effectively, the isolated imper- 
meable membrane approximation (m = 0) for A(q). 

Finally, Ramaswamy predicted that a dilute Helfrich- 
stabilized lamellar phase collapses when subjected to a 
flow field v — jyx in the a orientation jlTj . He considered 
the motion of a membrane in shear, and argued that only 
the confined fluid relaxation mode (m = 2) is relevant for 
wavelengths less than L p ||| : 

d t h + jyd x h = A V 2 ^+£(x,t). (2.12) 

The linearity of the advection term leads to an analytic 
form for the height correlation that decreases with shear. 
Ramaswamy demonstrated that at a critical shear rate 



(2.13) 



the fluctuations are significantly suppressed, provoking 
a layer collapse. Experimental confirmation of this was 
subsequently reported by Alkahwaji and Kellay [pTsj . 
Many other studies have incorporated the simple advec- 
tion term in the dynamics for u: for example, Cates and 
Milner studied the effect of flow on the isotropic-lamcllar 
transition |ll| , and Fredrickson studied the effect of flow 
on diblock lamellar phases [^0) . All of these treatments 
are suitable if flow does not significantly perturb the layer 
microstructure. In contrast the perturbation of the mi- 
crostructure in thermotropic lamellar phases was studied 
in Ref. || . In this case, shear flow was shown to introduce 
a tilt in the layers of a smectic-A liquid crystal. This led 
to layer reorientation and the possibility of an instability. 



Membrane in Flow in the c orientation; non-linear 
advection 



Now we consider the effect of a shear field v = 'yzx on 
a fluctuating membrane in the lamellar phase in the c ori- 
entation (Fig. |J). If we affinely transform the membrane 
in a small time 8t, the height field advects according to 
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• A 

FIG. 3: A fluctuating membrane at rest (solid line) and subject to a shear field (dashed line). 



h(x, y,t + St) = h(x — jzSt, y, t) 



^5t = ^(-jzdt) + 0(St) 2 . (2.14) 



The key point is that, because shear is assumed to advect the membrane, z should be equal to the height field h. This 
leads to a non-linear advective term, so that the equation of motion for a membrane in (d — 1) dimensions becomes: 



d t + A fh (x,t) 



d_ 

dx 



h(x,t) 



x'A (x - x' 



SH 



+ f(sc,t) 



-iojh(q, oj) + 17^] (g x — k x )h(k, tt)h(q — k, ui — O) 

I 



Sh(x',t) 
A(qWh(q)+aq,u). 



n k 



(2.15a) 
(2.15b) 



The advective term stretches the membrane, leading to 
a restoring "tension". In contrast, the relevant advec- 
tive term for a thermotropic in the c orientation re- 
mains ydu/dx, because the perturbation of the mem- 
brane structure is negligible, and undulations of u are 
always presumed to be of much smaller wavelength than 
the layer spacing. The noise has variance 

(£,(x 1 ,ti)£(x2,t 2 )) = D(xi - x 2 )S (xi - x 2 )5 (h -t 2 ) 

(2.16a) 

(£(<7i,^i)£(<Z2,w 2 )) = D(q 1 )6 d - 1 (q 1 + q 2 )6(u> 1 +w 2 ) 

(2.16b) 



D{q) 



2TA(q) 
Do\q\ m , 



2TA \q\ 



(2.16c) 
(2.16d) 



where Eq. ( 2.16c ) is the Fluctuation Dissipation Theo- 
rem. 

A similar equation was derived for a diffuse interface 
between phase separated domains subject to shear by 
Bray and co-workers |l(], O. However, there are two 
important differences. First, the "interface" of a mem- 
brane in the lamellar phase is intrinsically sharp. Second, 
although both bending (~ q 4 ) and surface (~ q 2 ) energies 



are present for a diffuse interface, at long length scales 
only the latter is relevant. On the other hand, surface 
tension is not present a priori in equilibrium lamellar 
phases, so the Helfrich Hamiltonian is dominated by the 

bending energy. 

For no flow 7 = 0, the non- linear term in Eq. ( 2.15b| ) 
vanishes and 



h(q,uj) = G(q,u)£(q,u) 



(2.17) 



This defines the bare linear propagator, equivalent to the 
linearized equation of motion, 



G- 1 (q,u) = -^ + A(q)S- 1 (q), 



(2.18) 



where S~ 1 (q) — ft|<j»| 4 . The poles of G(q,uj) yield the 
dispersion relation for the decay times of height fluctua- 
tions. The equal time height correlation function may be 
calculated from 



C(q) 



, 2n 
° dw 
2n 
= TS(q). 



(h*(q,u>)h(q,u>)) 



(2.19a) 

(2.19b) 
(2.19c) 
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Equations ( |2.15b| , |2 16bj ^.16d| ) define our model, with 



A(q) given by Eq. (£ 



In the next section we exam- 
ine the effect of the non-linear advective coupling on the 
spectrum of height fluctuations, through its effect on the 
linear propagator. 



III. COARSE-GRAINING PROCEDURE 

A. Description of the problem 

We wish to calculate the effective d ynamic al response 
of the height field, determined by Eq. ( |2.15b| ), for differ- 
ent models specified by the relaxation function A(q) and 
hence m. In Section [II B below we show that a one-loop 



perturbation analysis leads to a divergent response for all 
candidate values for m. In order to avoid this unphysi- 
cal divergence we ap ply a Renormalization Group (RG) 
procedure in Section III C to coarse-grain the system by 
removing the small scale and faster degrees of freedom, 
leaving an effective long wavelength theory. This proce- 
dure naturally generates a dynamic response that scales 
as q 2 , due to the advective non-linearity, is suggestive of 
a tension, and restores non-singular behaviour. 

In the study of dynamical critical phenomena |2l| the 
description of a system is coarse-grained until fixed points 
are found for which the system is self similar. Hence, 
a condition is found on the parameters of the dynami- 
cal equations of motion that leaves the dynamics of the 
height field, scale invariant. Contrary to a system at 
its critical point, the lamellar system here can only be 
coarse-grained up to the natural physical cutoff of the col- 
lision length. However, since this length is not necessarily 
in the scaling regime, such a calculation requires detailed 
knowledge of the flow equations for all new terms gener- 
ated in the equation of motion, which is generally impos- 



sible. On the other hand, if the collision length is close 
to, or larger than, the wavelength at which the scaling 
regime applies, we may use the simpler fixed point calcu- 
lation as a guideline to estimate the effective long wave- 
length dynamics. In fact, even if we are not in the scaling 
regime, the other terms generated by coarse-graining are 
of higher order in wave vector and are irrelevant at the 
fixed point, and are thus in the process of being driven 
to zero during the coarse-graining procedure. In the hy- 
drodynamic limit such terms are neglected. 

In the following analysis, we coarse-grain the system 
to generate the lowest order (in wave vector) additional 
term in the equation of motion, which depends on the 
coarse-graining length L p and the particular relaxation 
mechanism A(q), parametrized by m. However, unless 
the appropriate range for the length scales of a given re- 
laxation mechanism (see Fig. |J) encompasses the entire 
coarse-graining range, a full calculation demands a more 
precise choice of A(q). We shall not attempt such a cal- 
culation, but give the results under the assumption that 
a single value of m applies throughout the wave vector 
regime of interest. 

A consequence of the particular energy (~ q 2 ) intrin- 
sic to a diffuse interface is that Bray and co-workers were 
able to extract results only for "models" m > 0. In com- 
parison, for a fluctuating membrane (~ q 4 ), we may ex- 
amine the relaxation mechanisms m > — 2 including the 
case of an isolated impermeable membrane that relaxes 
by the hydrodynamic interaction of the surrounding sol- 
vent. 



B. One-loop correction to G(q,oj) 



We start by rewriting Eq. (2.15b) as 



h{q,uj) = G(q,uj) 



£(q,u) - ii^^iVx - k x )h(k,n)h(q -k,0J-tt) 



n fc 



(3.1) 



Eq. ( |3.l[ ) is shown in Fig. |J(a). Adding a perturbation f(q,oj) to Eq. (3^1) and averaging over the stochastic noise £ 
results in a renormalized linear propagator Gr, defined by 

d(h(q,cj)) 



G R (q,uj) = lim 

f^o df(q,uj) 



given by the Dyson equation (Fig. [I]b) p2] |, 

G R (q, us)- 1 



G(q,^)- 1 ~E( q , W ). 



(3.2) 



(3.3) 



The self energy E(q,w) shifts the poles of the effective propagator Gr and equivalently renormalizes the effective 
equation of motion (Eq. 2.18) for the long wavelength degrees of freedom. The one- loop contribution to the self 
energy E(q,w) (Fig. He) yields 



Q k 



X 



(q x - k x )G(-k, ~n)D{k) + k x G(k ~q,Vt- uj)D(k - q) 



(3.4) 



(a) 



h(q,a>) 

> 



(b) 



G(q,co) £(q,CO) G(q,co) 

— > X + — 



< 



iy EE(q-ig 



£2 k 




G R (q,co) 



>- 

G(q.fo) 



G(q,co) E(q,co) G/q,co) 



(c) 



q/2+k,co-£2 



-(q/2+k)-(co-£2) 



q,co 




q,co 



q,co 



q/2-k,Q 



q,co 



FIG. 4: (a) Diagrammatic representation of Eq. (3.1). Solid circles • represent non- linear vertices, (b) Dyson Equation for 
the renormalized propagator in terms of the bare propagator (single lines) and the self-energy (shaded circle), (c) One- loop 
correction to the self-energy. Open circles o represent noise contractions, (££}. 



A change of co-ordinates k — > q/2 + k and k — > q/2 — k in the first and second terms respectively of Eq. (3.4) gives 



n k 



G(% 



G[^ + k,uj-Q, 



Do 



shown in Fig. On symmetry grounds, the self energy can be written as 

S(q,w) = a 2 (q,uj)ql + a 4 (q,u)q* + 
in which case the renormalized propagator is given by 

Gr^uj)- 1 ~ + A(q)S'- 1 (q) +a 2 (q,w)^. 
and the associated noise correlation is 

(f (qi, (ga, wa)> = (D(gi) + D x q 2 x )S{ qi + q 2 )5{uj 1 + u 2 ). 



(3.5) 

(3.6) 
(3.7) 
(3.8) 



The effect of fluctuations on the slow hydrodynamic regime is due to the leading long wavelength static behavior, 
u> = 0. Hence, we expand to order q x for ui — > 0. Note that we implicitly ignore renormalizations of the frequency 
dependence, which should be sufficient to obtain scaling relations. 

After converting the sum to a continuum integral and integrating out the frequency (Eq. B5 in Appendix |b|) , the 
lowest order term in q x that appears in the static self-energy is 



£(g,0) 



• 2 2 

-7 q x 



D 



2(kA ) 



50 k d _^dkSd-i 
(2iT) d - 1 



1 



4|fcl 8 +' 



1k x (k x + k x kj_) 
2|fc| 12 + m 



(3.9) 



where S^-i is the unit sphere surface area in (d — 1) di- of our relaxation regimes m = —1,0,2, this integral di- 
mensions. For 10 + m — d > 0, which encompasses all 
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verges. In fact, before changing the lim its of integration 
in Eq. (3.4), the divergence in Eq. (3.4) can be traced to 
the second term in square brackets at q = k. That is, 
the divergence is actually independent of the lower cut- 
off. Such a singular response suggests that a description 
in terms of the microscopic degrees of freedom is physi- 
cally inconsistent. To proceed, we note that a tension-like 
term A x gJ added to the propagator yields non-singular 
behavior (Eq. |B7| in Appendix |b|). Furthermore, we will 
show below that the gradual thinning of degrees of free- 
dom during coarse-graining directly gen erates terms that 
render the response physical (Eq. B8 in Appendix ^). 
Hence, we argue that the correct physical description of 
the highly fluctuating system in shear must formally be 
derived by projecting out the small scale degrees of free- 
dom. We will outline this procedure next. 



C. Renormalization Group (RG) analysis 



ically, the original cutoff is A = n/a, where a is of or- 
der a molecular (surfactant) size, and we are interested 
in coarse-graining from this length to some, unknown, 
length, at which non-trivial scaling behavior is seen. The 
steps are as follows: 

1. The first step is to divide the Brillouin zone k G 
[0, A] into two parts: high wave vectors fc > E 
[A/6, A] to be removed, and the remaining long 
wavelengths fc < E [0, A/6]. The elimination of (as- 
sumed) fast modes results in an effective renormal- 
ized propagator G^(q,u>). Since there are no sin- 
gularities in this range of integration, only finite 
corrections to the parameters result. 

2. After coarse- graining, the resulting equation has a 
cutoff A/6. This difference from the original model 
is removed by rescaling the length scales, x, x±, h 
and the time scale t. 



Our goal is to successively "integrate out" the small 
scale fast degrees of freedom from the equation of motion, 
to yield an effective equation of motion for the remaining 
long wavelength degrees of freedom |23fl . Schematically, 
we can write 

h(x,t) = h>(x,t) + h < (x,t), (3.10) 

where h < (x,t) and /i > (x,t) are, respectively, small and 
large wave vector degrees of freedom. Upon removing 
the faster /i > (a;,t), an effective equation of motion for 
ft < (a;, t) will be generated. The form of the equation will 
differ from the original equation because of the non-linear 
advective term, which couples different modes together. 
Note that, generally, there are also non- linear terms (of 
order h 3 ) due to deviations from the Monge gauge limit 
and higher order advection terms, but the restriction of 
transverse length scales to within a patch length ensures 
that we remain close to this limit. The new equation of 
motion is most easily cast in terms of the renormalized 
propagator G R of the long wavelength degrees of free- 
dom. As noted above, we expect the contributions to the 
propagator, and generally to the noise, to be even pow- 
ers in q%, because of the symmetry of the non-linearity. 
In the hydrodynamic and long wavelength limit, and in- 
deed because higher order terms are irrelevant at the non- 
trivial fixed point, we focus on corrections of order q x . 

Essentially there are three steps to a momentum-shell 
RG in which, for convenience, we impose a short wave- 
length ultraviolet cutoff A in the ir-direction only; these 
are the fluctuations directly suppressed by flow. Phys- 



3. Finally we look for the fixed points of the recursion 
relation at which the theory is invariant under the 
first two steps. 

This procedure generates a recursion equation that 
may be used to find the behavior of the system in a scal- 
ing regime. Generally, this scaling regime corresponds to 
a description of the system at wavelengths longer than 
that wavelength at which the dynamics has effectively 
been driven to the fixed point. We will use this fixed 
point as an estimate of the dynamics of the system at 
wavelengths larger than the collision length L p . 

We follow Bray et ai. 's study Jl0|, [ll[] of the influence 
of shear flow on interfacial dynamics in a phase separat- 
ing system, which is governed by square gradient terms 
rather than quartic energy possessed by membranes. The 
scale transformation takes the form 

x = bx, x ± =b c x' ± , h = b x ti, t = b z t'. 

. (3-11) 

Since shear suppresses the fluctuations in the x-direction 
we expect to find £ < 1 when the shear is relevant. We 
will see that condition is only satisfied if m > — 2. Since 
we consider only models m > — 1, this condition is always 
satisfied. In such cases the transverse part q± dominates 
q x in the terms involving powers of \q\ so that the bare 
propagator is renormalized to 

G R l (q,u) = -iu + k{ qi _)S- l ( qi _)+K x ql (3.12) 

and the noise correlator is 



(£(<Zi, u; 2 )) = (D(q lx ) + D x q 2 x )5{q 1 + q 2 )S(cj 1 + uj 2 ). (3.13) 

I 

We have included the lowest order correction to Gr and the noise from the non-linearity. Applying the rescaling 
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Eq. (3.11) yields rescaled parameters in the equation of 
motion and the noise correlator: 



A' = b z -^ +m ^A 
A' x = b z - 2 A x + ... 

D' = b z-2x-l-?nC-(d-2) Do 

D' = b z - 2 x-^ d - 2 ^D x + . 



(3.14a) 
(3.14b) 
(3.14c) 
(3.14d) 
(3.14e) 



The parameters A^, and D' x acquire perturbative correc- 
tions due to the coarse-graining step of the RG proce- 
dure. In contrast 7, A and Dq do not acquire perturba- 
tive corrections. The nonrenor maliz ability of 7 follows 
from Galilean invariance of Eq. ( 2.14 ), which transforms 
t — > t + St and x — > x — -jhSt in the equation of motion. 

We first examine the linear theory to identify the crit- 
ical dimension d c . Since there are no perturbative cor- 
rections to a linear theory, t he requireme nt that the ex- 
ponents for b vanish in Eqs. (3.14b-3.14e) yields the con- 



ditions 



z = 2, 



Co 



Xo = 



\-m-2d 
2(m + 4) 



(3.15) 



The subsc ripts d enote the application to the linear the- 
ory. Eq. ( 3.14a ) determines the relevance of the shear 
rate 7 on t he co arse-graining, at the trivial fixed point. 
From Eq. (3.15) we obtain xo + z a — 1 = (jn — 2d + 
16)/ [2(4 + to)]. Therefore, 7 is relevant for d < d c where 



d r . = 



16 + to 



to > -2. 



(3.16) 



We can coarse-grain the theory perturbatively in 
Fourier space near the critical dimension d c of the the- 
ory. For d < d c we expect a new fixed poi nt to appear 
at which 7, A and D Q are non-zero. Eqs. ( 13.144 |3.14h| , 
3.14d) give the corresponding exponents exactly, 



J 



3(4 + to) 
14 + 2to - , 



to — 2d 
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From Eqs. ( 3.15 ) and ( 3.17 ) we see that £ < 1 for to > 
—2, in which cas e the approximation |q| ~ |qfj_| is consis- 
tent. From Eq. ( |Tf ) we find D' x = fe-8/(4+"0 D x , indi- 
cating that D x flows to zero at the fixed point (Fig. ||). 




Trivial 
Fixed Point 



Fixed Point 

X 



FIG. 5: Schematic flows of A. x and D x as a function of the 
coarse- graining length A. 

To find the values of A x and A at the non-trivial fixed 
poin t we now return to the RG procedure. Integrating 
Eq. (3.3) over the short wavelength modes gives the equa- 
tion for the effective renormalized propagator 



(3.18) 



2m — d ' 
I 



2(4 + to) 



(3.17) 



where E(q, u>) is given in E q. ( 3.5). Th en, se tting b = e l 
with / infinitesimal, Eqs. (3.14c) and (3.18) generate a 
differential flow equation for A x 



dA x 

~~dT 



A, 



(z-2) 



lim - — ^-r 



E(g,0) 



(3.19) 



Step (3) of the RG analysis consists of finding the fixed 
points A* for which the theory is invariant. This proce- 
dure is carried out in Appendix [b| for each value of to. 
Finally, we can make contact with our original discussion 
of an induced "tension" a, and extract a tension accord- 
ing to Eq. ([lj]) • The results are collected in Table Q. For 
771 = and to = 2 the procedure is straightforward and 
the results (and A x ) are independent of the low-A; cutoff 
of the theory. However, for m = — 1 we must cut off the 
theory at k — n/L p , and hence we find a result for A^ 
that depends on L p . 

As noted in the introduction, on the basis of estimating 
the height correlation function from the first term in the 
perturbation expansion in Eq. (|3.l|) , we would naively ex- 
pect the "tension" to scale quadratically with the strain 
rate; this scaling was also captured by considering the en- 
ergetic cost of bending and stretching a single membrane 
||. However, an anomalous scaling A x ~ 7 £m ,e m ^ 2 is 
generated in the scaling regime for to = and to = 2. 
The case of to = —1 is different. Owing to the divergence 
in the perpendicular direction at the low cutoff due to the 
more violent fluctuations at long wavelengths, the scaling 
does not follow the expected pattern (e_i = 2/3) of the 
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other mechanisms. However, whether accidentally or not, 
the first term in its power series satisfies £_i = 2 so at 
low strain rates it cannot be considered to exhibit anoma- 
lous scaling; in general though, contributions from higher 
order terms indicate that the scaling is also anomalous. 

In addition, equating bending and tension energies 
leads to an expression for the shear rate j c at which fluc- 
tuations are significantly suppressed, 



r (m+3)/2 Ao 
K (m+l)/2gm+i ' 



(3.20) 



where in the case m = — 1 there are higher order con- 
tributions, shown in Table |. If we compare the critical 
shear rate for isolated impermeable membrane relaxation 
with Bruinsma's result (Eq. 2.11), the scaling is the same 
apart from a factor of yT '/ ' k. 



IV. RESULTS AND DISCUSSION 

A. Single Membrane Dynamics in the Scaling 
Regime 

We have considered the dynamics of a single mem- 
brane in the c-orientation (Fig. |l|) with respect to shear 
flow. Advection couples different Fourier modes, and 
hence renormalizes the effective response. We have es- 
timated the effective long-wavelength theory, that would 
be obtained by removing smaller and faster degrees of 
freedom with wave vectors q > A, by calculating the be- 
havior of the fluctuating membrane in the scaling regime. 
This dynamic coarse-graining generates, to lowest order, 
a term —A x q 2 in the long wavelength propagator. This 
is the principle qualitative result of this work. 

The function A^ depends on wave vector q, the qui- 
escent relaxation mechanism A(q), the strain rate, and 
the wavector scale A to which coarse-graining has been 
performed (Appendix 0) : 



D. One-Step Coarse-Graining 



A_ T ~ A™ A 



'7 



(4.1) 



So far we have demonstrated how to calculate the non- 
trivial scaling behavior of the membrane, assuming that 
fluctuations generate a tension-like term that renders 
contributions to the self-energy non-singular. Here, we 
show explicitly how the "first step" of a coarse-graining 
procedure produces such a term. Here we coarse-grain 
the system in "one step" , by making a small perturba- 
tion to the original microscopic cutoff A = n/a, where a 
is a typical molecular dimension. For a small perturba- 
tion about the trivi al fix ed point (A^ = 0) the differential 
flow equation Eq. ( 3.19 ) becomes 



dA x 



lim -^S(g,0). 



(3.21) 



We show at the end of Appendix [b] the calculation of the 
coarse-grained self-energy in the limit A x 
Inserting the resulting self-energy (Eq 



(Eq.|B20[). 
B2~l"l ) into the 
recursion relation yields an expression for A^ and, via 
Eq. (1.3), a "tension" that depends on the cutoff and the 
relaxation mechanism, 



« 2 AgA 9 + r 



- d rk~ 



(3.22) 



where f3 is a numerical prefactor and I is a small number 
that depends on the chosen coarse-graining step. As ex- 
pected, for all relaxation mechanisms, a first step coarse- 
graining leads to a "tension" that is the same as the 
naive scaling ~ i] 2 ^ 2 . Thus we may infer that the coarse- 
graining process modifies the dependence of the scaling 
for m = 0, 77i — 2 and also for m = — 1 (but not for very 
small strain rates). 



where m parametrizes the relaxation mechanism A(q) = 
Ao9 m , and the constant A m depends on T, k, L p , and Aq. 
The results are summarized in Table [j]. This restoring 
term is suggestive of an anisotropic "tension" a, which 
would appear in the dynamics as aq 2 h < {q) (Eq. |1.3[ ), 
with cr = A x /A(q), except that the non-analytic form 
generally leads to a wave vector dependence. In the per- 
meable limit the wavector dependence is absent, v = 0, 
while in other cases there is a wave vector dependence. 
Hence, referring to the newly generated term as a tension 
is suggestive at best. Nonetheless, this term can be ex- 
pected to suppress fluctuations, and hence influence the 
effective collision rate, and in turn the Helfrich interac- 
tion potential, in the presence of shear flow. Elsewhere, 
we have used an effective energetic tension to parametrize 
the reduction in fluctuations and the corresponding flow- 
induced strain or change in layer spacing j7|. 

It is important to recognize that, although the "ten- 
sion" in Table | applies, strictly, only to wavelengths of 
order the collision length, it is generated at all wave- 
lengths larger than the smallest cutoff and grows during 
the coarse-graining procedure. In reference j?j we re- 
placed this wavevector-dependent tension by an average 
value that applies for all wave vectors. This certainly 
changes any quantitative predictions, but does not influ- 
ence the qualitative aspects of those results. This naive 
estimate should evidently be replaced by a much more 
sophisticated dynamic analysis that simultaneously per- 
forms the dynamic coarse-graining in the presence of the 
advective non-linearity and a self-consistent (or coarse- 
graining) procedure to recover the Helfrich interaction 
behavior that stabilizes the lamellar stack. Such a calcu- 
lation is beyond the scope of this work. 

If the layer spacing does adjust in flow due to an in- 
duced tension, a non-Newtonian response is likely to be 
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TABLE I: Row (2) shows the fixed point s co rresponding to different membrane relaxation mechanisms, given in row (1). Row 
(3) gives the "tension", defined by Eq. jl.3j ) as a function of the cutoff, A and the relaxation mechanism. Rows (4) and (5) 
display the results for A -1 ~ L v and q _1 ~ L p ~ lyj k/T respectively. Row (6) displays the critical shear rate for suppression 
of fluctuations. 
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found. Most probably this will be shear thinning, be- 
cause of the greater local regularity of the flow, although 
it is not obvious that this is the case. The magnitude of 
the viscous response is a complicated balance of dissipa- 
tion incurred within bilayers, and local inhomogeneous 
shears due to the fluctuating layers. The single study 
that reported a change in layer spacing also reported a 
shear thinning response Q. Shear thinning behaviour 
has been observed in some Helfrich-stabilized systems in- 



cluding C 12 E 5 §, AOT 1 24 and SDS p§. 



B. Effective Long Wavelength Dynamics 

The effective long wavelength dynamics of the single 
membrane is of the form (in Fourier space) 



-iuh<(q, u) + ~ k x )h<{k, n)h<(q - fe, w - fi) = - [HqW + A x q 2 x ] h<(q) + £<(q, w), (4.2a) 

n k 

(qi, wi) (92, wa)) = (-D(gi) + D x ql)6{ qi + q 2 )6(uj 1 + u> 2 ), (4.2b) 

where /^(q, w) is the small wave vector (coarse-grained) height field, and the noise in principle, incorporates the 
eliminated degrees of freedom in addition to the original small scale degrees of freedom. This yields a proportionality 
between correlation and response, and a generalized fluctuation-dissipation theorem (FDT) is satisfied, although the 
simple proportionality factor of temperature relating correlation and response is replaced by the more complicated 
noise correlations. In the case, which we have assumed, that the scaling limit is reached before the patch size has 
been reached, D x vanishes and an effective temperature, albeit shear rate dependent, can be ascribed to the system 
according to the fixed point value for D. 

Ideally, coarse-graining should continue until all wavelengths less than the collision length L p have been removed, at 
which point the resulting theory would be used as a starting point for understanding the dynamics of the usual mean 
smectic layer displacement u, rather than the microscopic membrane position h. Note that, at this point, collisions 
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intervene in a non-trivial way to limit affine layer advection, and the coarse-grained smectic phase variable u advects 
according to jyd x u rather than jud x u. This behavior should, in principle, emerge smoothly in an ideal calculation. 

The resulting dynamics of a strongly fluctuating layered system in shear flow are best cast in terms of the velocity 
field, in the standard two-fluid form [g6|, as 



p (d t + v ■ V) v = -Vp + ?/V 2 t> + n f n 
(d t +v ■ V)u = v z , 



(4.3) 
(4.4) 



where we have, for convenience, shown the form in the absence of permeation. The normal force f n differs from the 
usual normal force by the term generated upon coarse-graining, 



fn 



6u(q,t) d 



(4.5) 



where the free energy J- should also include the layer 
compression energy density ^B(d z u) 2 . Note the factor 
of d in the second term, reflecting the inherent three di- 
mensional nature of smectic elasticity. The noise defines 
an effective temperature that is generally not the phys- 
ical temperature, and may have additional correlations 
D x q x that reflect the flow (depending on whether or not 
the scaling regime has been reached). 

The additional term is only present for strong flows, 
and penalises layer undulations in the x direction; this 
is because such undulations are performed at the ex- 
pense of the microscopic height fluctuations, which are 
highly stretched in strong flows. This term is not ex- 
pected to appear in situations where the microstructure 
of the smectic layers is essentially undisturbed by flow, 
as in typical thermotropic smectics (but see the calcula- 
tion of Auernhammer et al. p] for a counter example). 
One could also envision this term as a non-equilibrium 
contribution to an effective free energy, which has been 
postulated by Jou and co-workers in their studies of com- 
plex fluids using extended irreversible thermodynamics 
J27J; however the dependence on strain rate that we de- 
rive, ~ 7 em , is not necessarily analytic, unlike their 
assumptions. 

It is important to remember that the generation of 
the tension-restoring term is only one of several possible 
dynamic effects; other effects include the rearrangement 
of defect distributions, which is also likely to lead to a 
shear thinning response 0] . 



C. Summary 

In this work we have studied the effect of flow on the 
dynamics of fluctuating membranes. We have made sev- 
eral assumptions, which we collect for completeness: 

• We assumed that an £-expansion is sufficient to 
describe the effect of coarse-graining the theory 
up to the collision length; in this limit the renor- 
malized noise reduces to an effective temperature. 
Whether or not scaling is truly reached is an open 



question. It is more likely that there are resid- 
ual noise correlations when the collision length has 
been reached. Moreover, the critical dimension d c 
is quite high and fluctuations are quite important; 
we have considered m = — 1, 0, 2, for which, respec- 
tively, dc = 8, 9. 

• Since an e-expansion is not likely to hold so far from 
the critical dimension, our calculation is strictly a 
self-consistent one-loop calculation. 

• We have considered the different membrane relax- 
ation mechanisms (permeable, squeezing, isolated) 
separately. In reality, the mechanism changes dur- 
ing the coarse-graining process, according to Fig. ^; 
nonetheless, this does not detract from our primary 
message, that the perturbation of the microstruc- 
ture of highly fluctuating membranes can lead to 
an additional restoring term in the long wavelength 
dynamics. 

• The coarse-graining can be performed only up to 
L p , because at this length scale the long range steric 
repulsion is important. In fact, we have completely 
ignored steric interactions. A more precise treat- 
ment would involve simultaneously treating flow 
and collisions, or treating the flow within a self- 
consistent scheme using, for example, a harmonic 
potential to mimic collisions. 

The significant accomplishments of this study have 
been, first, a qualitative estimate of the effect of flow 
on highly fluctuating lamellar phases. More importantly, 
however, is the demonstration that flow can strongly 
modify the fluctuation spectrum and generate new ef- 
fects in the macroscopic response, via an RG-like self- 
consistent coarse-graining technique. The theory that 
emerges has a natural effective noise that need not sat- 
isfy the usual equilibrium fluctuation-dissipation theo- 
rem. Similar renormalizations of hydrodynamic descrip- 
tions can be expected for other complex fluid systems 
with highly fluctuating mesoscopic degrees of freedom, 
such as wormlike micellar systems (for example, the mi- 
cellar length, modulus, and relaxation times could be ex- 
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pected to renormalize due to the effect of flow on the 
local charge distribution and undulation spectra). 
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APPENDIX A: PERMEATION LENGTH SCALE 

We estimate the permeation length £ by assuming that 
permeation is dominated by solvent flow through pores 
in the membranes; indeed, at certain surfactant concen- 
trations pores are very common - see for example the 
studies [|9[ pQ] . We assume cylindrical pores of diame- 
ter 2w within a membrane of thickness r, separated by a 
mean distance R (Fig. ^|). We wish to derive the kinetic 
coefficient Aq for layer relaxation, defined by 



dth 



-A 



5H 
] ~Sh' 



(Al) 



Identifying the mean solvent flow velocity (v) with the 
relaxation speed of the membrane dth yields 



(v) ~ A 



5H 
~6h' 



(A2) 



The pressure difference Ap = p + —p~ driving solvent flow 
is given by the force per area STi/Sh that the membrane 
exerts on the fluid. Thus the kinetic coefficient is 



A ~ {V) 
Ao A~p 



(A3) 



For each pore the mean velocity is the flux of material Q 
flowing through a pore per unit surface area of membrane 
A, 



(v) = 



Q 

A' 



(A4) 



We note that in terms of permeability V, the flux is given 
by Q/A = -Wp/r) (Darcy's law). Thus V ~ A ?7T, 
which was used by Leng |3lj in the context of swelling 
compressed lamellar phases. To calculate the flux we 
assume that the pressure gradient sets up a Poiseuille 
flow given by the viscous flow force balance 



Ap 

T 



(A5) 



where v c is the velocity of the solvent at the centre of the 
flow, and so 



Ap 



A R 2 t Tj 



(A6) 



On s ubs tituting the express ion for the mean velocity in 
Eq. ( |A3| ) by the flux in Eq. (At) we obtain an expression 
for the kinetic coefficient 



A 



w 



-R 2 ri' 



(A7) 



Therefore, the "permeation length" is C — w i /(R 2 r). 
APPENDIX B: CALCULATION OF THE 
SELF-ENERGY AND THE CRITICAL POINTS 
FOR DIFFERENT RELAXATION MECHANISMS 



In this appendix we calculate the coarse-grained self-energy S(q, 0), which we insert into the recursion relation 
(Eq. 3.19) to find the fixed points A* for different values of m. We also calculate the general expression for the 
self-energy, and the associated A x derived from the recursion relation in a "first step" coarse-graining process. 



1. Coarse-grained self energy 



After demonstrating that a naive calculation of the self-energy, equivalent to A^ = 0, leads to a divergence, the 
main task here is to show ho w the divergence is eliminated for A x ^ 0, enabling us to proceed with the integration. 
Having established in Section IIIB that the self energy (Eq. 3.5) can be written in powers of q x , we commence from 
its expression in terms of the renormalized propagator (Eq. 3.7) and its corresponding noise correlation (Eq. tL 



E(q, u) = -7 2 ££^!- fc *) G *(f- fe > fi 



Gr 



k,u> — n 



At 



k 



2 



(Bl) 



We consider the slow hydrodynamic regime u> — > 0, and ignore frequency-dependent corrections, so that 

G R (q, uj)- 1 ~ -iw + A(q)S- 1 (q) + A x q 2 x , 

where A x = 02(9,0), A(q) = Aoq m and S^ 1 (q) — nq 4 . The disregard for frequency-dependent corrections should 
suffice to obtain scaling properties. The sums are converted to integrals by Eq. (2.2) 

dn d^kSa-i 
"2T 



EE 



■n/L p 



(27T) 



d-1 



(B2) 



(B3) 
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FIG. 6: Pore defects in a membrane (side view). 



where S n = n-n n l 2 /T(n/2 + 1) is the surface area of an n dimensional unit sphere and the limits in k are given by 
the physical cutoffs. The first step is to perform the fi-integral by contour integration, yielding a positive pole at 

in = - \a x (k x + 3f) 2 + kA q \k + f | 4+m | , 



E(g,0) = -y 



k d - 2 dkS d -i 




D \k + % | m 4 


- Dx (k x 


+ f) 2 




(2 7 r) d - 1 


K (k x + + kAq \k + | 4+m 




~2A X (kl + & 


+ kA 


fc+f| 


4+m + K A |fc-f| 4+m 



(B4) 



As we explained in Section [II B, the leading long wavelength behaviour arises from an expansion of Eq. (B4) to 
lowest order in q x , which we now investigate for the cases A^ = and A^ ^ 0. For A x — 0, and thus also D x — 0, the 
self-energy becomes 



S(q,0) = -7\ 2 



At 



2(kA ) 



50 k d - 2 dkS d -! 



:d-l 



1 



4|fcl 8 +' 



2|fe| 12 + m 



(B5) 



which diverges at low k for 10 + m — d > (note that all relaxation mechanisms we consider obey m > d — 10). As 
we discuss at the end of Section III B , this divergence is unphysical. Upon coarse-graining the theory, the term A x q x 
will be generated, which obviously changes the character of the integral. Next we show that the implementation of 
coarse-graining such that A x ^ removes this divergence. 



Following the first step of the renormalization group analysis in Section 1IIC the removal of high wave vectors fc > 
in the x direction in the range \e~ l < k x < X is equivalent to a change in limits, 



E 

fe> 



2tt 



oo jd-2 



k\S d 



(2tt 



)d-2 



(B6) 



Note that there is no change in the limits in the perpendicular direction. We showed in Section |III C| th at when 
A x 7^ and m > — 1, we approximate |fe| — > |fej_| so that the renormalized propagator becomes Eq. (3.12) and the 
noise becomes Eq. ( |3.13 ). In addition, within the scaling regime or close to the fixed point, D x — > 0. With these 
assumptions the resulting approximation to the lowest order expansion of the self energy is 



S(g,0) 



oo jd—3 



dk iS, 



±Jd-2 



k'. 



(2tt 



id-2 



{A x kl + KA k 



1+m) 



A x k x 



2 {A x kl + KA k 



4+m\ 



0(9 4 )- (B7) 



The divergence that we encountered before for A^ = is now eliminated, which enables us to proceed with the 
integration over all wave vector s in t he transverse direction, in order to calculate the self-energy E(q, 0) to be inserted 
into the recursion relation (Eq. |3 . 1 £ ) . 



2. Calculation of the fixed points 



With the condition that d + m— 2 > (a criterion that we discuss later), integration of Eq. (B7) leads to 

(d + m - 2 > 0) , (B8) 



\ (14 + 2m-cf) , 

S(q,0) = ~ 2 ^ W-d UAx<lJ 
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where 



U = Y 



S d - 



d + m- 2\ ^ /14 + 2m -d 
\ m + 4 



8(m + 4) (27r) d - 2 
r(u) is the Gamma function and we have used 

x z - x dx _ T(z)T(w) 



Dq(kA ) ™+ 4 A x 



QR(z) > 0,%) > 0). 



(B9) 



(BIO) 



(i 



The negative exponent of in Eq 
for coarse-graining. 



(l + x) z+w T(z + w) 

reveals the divergence established earlier and therefore reaffirms the necessity 



Inserting the self energy (Eq. jB8|) into Eq. ( p,19[) leads to an expression for the recursion relation for , which is 
most conveniently written in terms of the coupling constant U: 



dU 
~dT 



16 + m- 2d 



U 



(14 + 2m - d) 2 2 
(m + 4)(m + 10-d) ' 



(Bll) 



where we have used Eq. (3.17) to eliminate z. Having thus removed the high wave vectors and rescaled all the 
parameters the final step of the RG analysis is to find the fixed points of Eq. (Bll) for which the theory is invariant. 
Consistent with the previous determination of the critical dimension, the linear term changes sign for d = d c = 
(16 + m)/2. Since the quadratic term is negative for d < m + 10 (or d < d c for m < —4), there is a non-zero stable 
fixed point 



U* 



9(4 - 6) 



(B12) 



to first order in e where e = d — d c . This shows that the RG perturbation is well-behaved and the exponents are 
correct. For d > d c the only stable fixed point is U* — 0, corresponding to an irrelevant non-linearity and recovering 
the exponents for the linear theory. Thus in the general case for d + m — 2 > 0, 



where 



A! 



Otm+l 



a 



m+l 



I -m^.2— d-m^— 16— m+2d^-2m-8^2(4+m) 



T 



l/(14+2m-d) 



4(m + 4)(27r) 



d-l 



16 + to- 2d 



d + m-2\ (14 + 2m- d 
\ 4 + to 



(4+m)/(14+2ro-d) 



(B13) 



(B14) 



We may apply this result to two of the relaxation mechanisms that we considered; both results are found in Table || 
in Section IV. For the permeable case, d = 3 and m = 0, 



A* 



Cti 



rp4/U A 3/11 

1 A ^8/11 

3/11 ' ' 



K 1/11 A 10/ 

For the confined case, d = 3 and m = 2, 

T 2/5 A 4/5 



where 



ai 



189Si (1 
2560?r 2 V 4 



-i 4/11 



0.342 



(B15) 



A* = a 3 — 



K l/5 A 4/5 



^4/ 5) 



where 



a 3 



-i 2/5 



0.378 



(B16) 



256tt 2 

A similar analysis cannot be cond ucte d for the case d = 3, to = — 1. In the hydrodynamic limit a divergence in the 
lower limit in the k± integral of Eq. (B7) occurs for d + m — 2 < 0. By introducing a lower cutoff given by the inverse 
collision length of the system L~ 1 and writing the integral in terms of the dimensionless quantity y = L^A X X 2 / 'kAo, 
the expression for the self energy for d = 3, to = — 1 becomes 



E(q,0) 



ao7 2 A) 
A 3 A3 



3 Mn(y+1)- 



y + 1 



where 



a 



Si 
96tt 2 



0.000528. 



(B17) 



We now proceed as before by combining Eq. ( 3. IS ) and the expression for the self-energy (Eq. B17 ) to give the 
differential flow equation in terms of y, 
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ao7 2 -Po 
A 3 



^A 2 
nA y 



3 ln(y + 1) 



y 

(y + i) 2 



(B18) 
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The unstable fixed point corresponding to the irrelevant non-lin earit y is given by = 0. As we are unable to give 
an analytic expression for the stable fixed point solution of Eq. (BIS) we show instead its power series in j 2 , 



a: 



TL A . 2 

K 2 A 



(B19) 



Hence for small strain rates A* ~ 7 2 . 



3. First step coarse-graining 



Here we demonstrate the procedure for a "first step" coarse-graining of the self-energy in Section HID, i.e. we 
calculate A x by perturbing about the trivial fixed point A* = 0. First we return to the expression for the coarse- 
grained self-energy in Eq. (B4) derived from Eq. (Bl). Upon coarse-graining in the x direction, the limits of the 
sum and the integration are changed according to Eq. (B6). However, due to being far from the scaling regime, the 
assumption that \k\ — > |fej_| no longer holds. Hence in the limit A^ — > (and thus D x — > 0) appropriate for a "first 
step" coarse-grain, the integral becomes 



£( 9 ,0) 



• 2 2 

-7 Q x ; 
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dkS d -2 



(2n 



: 2(kA ) 2 

from which we may compute the self-energy for general m and d, 

T 



1 



S(q,0) = -/3- 



4|fe|8+m 



• 2 2 7 



2k x (k'^ + k x k^) 
2\k\ 12 + m 



V(qi 



K 2 A X 9+r 

(3 is a numerical prefactor that depends on the relaxation mechanism m and the dimension d, 



(B20) 
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(B22) 



On substituting Eq. ( B21| ) into the recursion relation for around the trivial fixed point (Eq. 3.21) we find that. 

T 



dk x {l)~p 



K 2 A X 9+m - d 



(B23) 



which leads to the "tension" that depends on the cutoff in Eq. (3.22). 
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